home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2001 May / SGI IRIX Base Documentation 2001 May.iso / usr / share / catman / p_man / cat3 / libblas / syr.z / syr
Encoding:
Text File  |  1998-10-30  |  13.5 KB  |  265 lines

  1.  
  2.  
  3.  
  4. ____SSSSYYYYRRRR,,,,____HHHHEEEERRRR((((3333FFFF))))                                                    ____SSSSYYYYRRRR,,,,____HHHHEEEERRRR((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      dsyr, ssyr, zher, cher - BLAS Level Two   (Symmetric/Hermitian)Matrix
  10.      Rank 1 Update
  11.  
  12.  
  13. FFFFOOOORRRRTTTTRRRRAAAANNNN 77777777 SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee ddddssssyyyyrrrr(((( uuuupppplllloooo,,,, nnnn,,,, aaaallllpppphhhhaaaa,,,, xxxx,,,, iiiinnnnccccxxxx,,,, aaaa,,,, llllddddaaaa ))))
  15.            cccchhhhaaaarrrraaaacccctttteeeerrrr****1111        uuuupppplllloooo
  16.            iiiinnnntttteeeeggggeeeerrrr            nnnn,,,, iiiinnnnccccxxxx,,,, llllddddaaaa
  17.            ddddoooouuuubbbblllleeee pppprrrreeeecccciiiissssiiiioooonnnn   aaaallllpppphhhhaaaa
  18.            ddddoooouuuubbbblllleeee pppprrrreeeecccciiiissssiiiioooonnnn   aaaa(((( llllddddaaaa,,,,****)))),,,, xxxx((((****))))
  19.  
  20.      ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee ssssssssyyyyrrrr(((( uuuupppplllloooo,,,, nnnn,,,, aaaallllpppphhhhaaaa,,,, xxxx,,,, iiiinnnnccccxxxx,,,, aaaa,,,, llllddddaaaa ))))
  21.            cccchhhhaaaarrrraaaacccctttteeeerrrr****1111        uuuupppplllloooo
  22.            iiiinnnntttteeeeggggeeeerrrr            nnnn,,,, iiiinnnnccccxxxx,,,, llllddddaaaa
  23.            rrrreeeeaaaallll               aaaallllpppphhhhaaaa
  24.            rrrreeeeaaaallll               aaaa(((( llllddddaaaa,,,,****)))),,,, xxxx((((****))))
  25.  
  26.      ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee zzzzhhhheeeerrrr(((( uuuupppplllloooo,,,, nnnn,,,, aaaallllpppphhhhaaaa,,,, xxxx,,,, iiiinnnnccccxxxx,,,, aaaa,,,, llllddddaaaa ))))
  27.            cccchhhhaaaarrrraaaacccctttteeeerrrr****1111        uuuupppplllloooo
  28.            iiiinnnntttteeeeggggeeeerrrr            nnnn,,,, iiiinnnnccccxxxx,,,, llllddddaaaa
  29.            ccccoooommmmpppplllleeeexxxx****11116666         aaaallllpppphhhhaaaa
  30.            ccccoooommmmpppplllleeeexxxx****11116666         aaaa(((( llllddddaaaa,,,,****)))),,,, xxxx((((****))))
  31.  
  32.      ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee cccchhhheeeerrrr(((( uuuupppplllloooo,,,, nnnn,,,, aaaallllpppphhhhaaaa,,,, xxxx,,,, iiiinnnnccccxxxx,,,, aaaa,,,, llllddddaaaa ))))
  33.            cccchhhhaaaarrrraaaacccctttteeeerrrr****1111        uuuupppplllloooo
  34.            iiiinnnntttteeeeggggeeeerrrr            nnnn,,,, iiiinnnnccccxxxx,,,, llllddddaaaa
  35.            ccccoooommmmpppplllleeeexxxx            aaaallllpppphhhhaaaa
  36.            ccccoooommmmpppplllleeeexxxx            aaaa(((( llllddddaaaa,,,,****)))),,,, xxxx((((****))))
  37.  
  38.  
  39.  
  40. CCCC SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  41.      vvvvooooiiiidddd ddddssssyyyyrrrr(((( uuuupppplllloooo,,,, nnnn,,,, aaaallllpppphhhhaaaa,,,, xxxx,,,, iiiinnnnccccxxxx,,,, aaaa,,,, llllddddaaaa ))))
  42.            MMMMaaaattttrrrriiiixxxxTTTTrrrriiiiaaaannnngggglllleeee     uuuupppplllloooo;;;;
  43.            IIIInnnntttteeeeggggeeeerrrr            nnnn,,,, iiiinnnnccccxxxx,,,, llllddddaaaa;;;;
  44.            ddddoooouuuubbbblllleeee             aaaallllpppphhhhaaaa;;;;
  45.            ddddoooouuuubbbblllleeee             ((((****aaaa))))[[[[llllddddaaaa****nnnn]]]],,,, ((((****xxxx))))[[[[ nnnn ]]]];;;;
  46.  
  47.      vvvvooooiiiidddd ssssssssyyyyrrrr(((( uuuupppplllloooo,,,, nnnn,,,, aaaallllpppphhhhaaaa,,,, xxxx,,,, iiiinnnnccccxxxx,,,, aaaa,,,, llllddddaaaa ))))
  48.            MMMMaaaattttrrrriiiixxxxTTTTrrrriiiiaaaannnngggglllleeee     uuuupppplllloooo;;;;
  49.            IIIInnnntttteeeeggggeeeerrrr            nnnn,,,, iiiinnnnccccxxxx,,,, llllddddaaaa;;;;
  50.            ffffllllooooaaaatttt              aaaallllpppphhhhaaaa;;;;
  51.            ffffllllooooaaaatttt              ((((****aaaa))))[[[[llllddddaaaa****nnnn]]]],,,, ((((****xxxx))))[[[[ nnnn ]]]];;;;
  52.  
  53.      vvvvooooiiiidddd zzzzhhhheeeerrrr(((( uuuupppplllloooo,,,, nnnn,,,, aaaallllpppphhhhaaaa,,,, xxxx,,,, iiiinnnnccccxxxx,,,, aaaa,,,, llllddddaaaa ))))
  54.            MMMMaaaattttrrrriiiixxxxTTTTrrrriiiiaaaannnngggglllleeee     uuuupppplllloooo;;;;
  55.            IIIInnnntttteeeeggggeeeerrrr            nnnn,,,, iiiinnnnccccxxxx,,,, llllddddaaaa;;;;
  56.            ZZZZoooommmmpppplllleeeexxxx            aaaallllpppphhhhaaaa;;;;
  57.            ZZZZoooommmmpppplllleeeexxxx            ((((****aaaa))))[[[[llllddddaaaa****nnnn]]]],,,, ((((****xxxx))))[[[[ nnnn ]]]];;;;
  58.  
  59.      vvvvooooiiiidddd cccchhhheeeerrrr(((( uuuupppplllloooo,,,, nnnn,,,, aaaallllpppphhhhaaaa,,,, xxxx,,,, iiiinnnnccccxxxx,,,, aaaa,,,, llllddddaaaa ))))
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. ____SSSSYYYYRRRR,,,,____HHHHEEEERRRR((((3333FFFF))))                                                    ____SSSSYYYYRRRR,,,,____HHHHEEEERRRR((((3333FFFF))))
  71.  
  72.  
  73.  
  74.            MMMMaaaattttrrrriiiixxxxTTTTrrrriiiiaaaannnngggglllleeee     uuuupppplllloooo;;;;
  75.            IIIInnnntttteeeeggggeeeerrrr            nnnn,,,, iiiinnnnccccxxxx,,,, llllddddaaaa;;;;
  76.            CCCCoooommmmpppplllleeeexxxx            aaaallllpppphhhhaaaa;;;;
  77.            CCCCoooommmmpppplllleeeexxxx            ((((****aaaa))))[[[[llllddddaaaa****nnnn]]]],,,, ((((****xxxx))))[[[[ nnnn ]]]]
  78.  
  79.  
  80.  
  81. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  82.      ddddssssyyyyrrrr and ssssssssyyyyrrrr perform the symmetric rank 1 operation
  83.  
  84.            A := alpha*x*x' + A,
  85.  
  86.      zzzzhhhheeeerrrr and cccchhhheeeerrrr perform the hermitian rank 1 operation
  87.  
  88.            A := alpha*x*conjg( x' ) + A,
  89.  
  90.      where alpha is a real/complex scalar, x is an n element vector and A is
  91.      an n by n symmetric/hermitian matrix.
  92.  
  93.  
  94. PPPPAAAARRRRAAAAMMMMEEEETTTTEEEERRRRSSSS
  95.      uuuupppplllloooo    On entry, uuuupppplllloooo specifies whether the upper or lower triangular
  96.              part of the array A is to be referenced a follows:
  97.  
  98.  
  99.  
  100.              FFFFOOOORRRRTTTTRRRRAAAANNNN
  101.              uplo = 'U' or 'u'   Only the upper triangular part of A
  102.                                  is to be referenced.
  103.              uplo = 'L' or 'l'   Only the lower triangular part of A
  104.                                  is to be referenced.
  105.  
  106.              CCCC
  107.              uplo = UpperTriangle     Only the lower triangular part of A
  108.                                       is to be referenced.
  109.              uplo = LowerTriangle     Only the lower triangular part of A
  110.                                       is to be referenced.
  111.  
  112.              Unchanged on exit.
  113.  
  114.      nnnn       On entry, nnnn specifies the the order of the matrix A.  nnnn must be
  115.              at least zero.
  116.              Unchanged on exit.
  117.  
  118.      aaaallllpppphhhhaaaa   specifies the scalar alpha.
  119.              Unchanged on exit.
  120.  
  121.      xxxx       Array of size at least ( 1 + ( n - 1 )*abs( incx ) ). Before
  122.              entry, the incremented array xxxx must contain the n element vector
  123.              x.
  124.              Unchanged on exit.
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.  
  134.  
  135.  
  136. ____SSSSYYYYRRRR,,,,____HHHHEEEERRRR((((3333FFFF))))                                                    ____SSSSYYYYRRRR,,,,____HHHHEEEERRRR((((3333FFFF))))
  137.  
  138.  
  139.  
  140.      iiiinnnnccccxxxx    On entry, iiiinnnnccccxxxx specifies the increment for the elements of xxxx.
  141.              iiiinnnnccccxxxx must not be zero.
  142.              Unchanged on exit.
  143.  
  144.      aaaa       An array containing the matrix A.
  145.  
  146.              FFFFOOOORRRRTTTTRRRRAAAANNNN
  147.              Array of dimension ( lda, n ).
  148.  
  149.              CCCC
  150.              A pointer to an array of size lda*n.
  151.              See note below about array storage convention for C.
  152.  
  153.              Before entry with uuuupppplllloooo ==== ''''UUUU'''' oooorrrr ''''uuuu'''' oooorrrr , the array elements
  154.              corresponding to the leading n by n upper triangular part of the
  155.              matrix A must contain the uppertriangular part of the
  156.              symmetric/hermitian matrix and the corresponding strictly lower
  157.              triangular part of A is not referenced. On exit, the array
  158.              elements corresponding to the upper triangular part of the matrix
  159.              A is overwritten by the upper triangular part of the updated
  160.              matrix.
  161.  
  162.              Before entry with uuuupppplllloooo ==== ''''LLLL'''' oooorrrr ''''llll'''' oooorrrr , the array elements
  163.              corresponding to the leading n by n lower triangular part of the
  164.              matrix A must contain the lower triangular part of the
  165.              symmetric/hermitian matrix and the corresponding strictly upper
  166.              triangular part of A is not referenced. On exit, the array
  167.              elements corresponding to the lower triangular part of the matrix
  168.              A is overwritten by the lower triangular part of the updated
  169.              matrix.
  170.  
  171.              Note that the imaginary parts of the diagonal elements need not
  172.              be set and are assumed to be zero.
  173.  
  174.      llllddddaaaa     On entry, llllddddaaaa specifies the first dimension of A as declared in
  175.              the calling (sub) program. llllddddaaaa must be at least ( k + 1 ).
  176.              Unchanged on exit.
  177.  
  178.  
  179. CCCC AAAARRRRRRRRAAAAYYYY SSSSTTTTOOOORRRRAAAAGGGGEEEE CCCCOOOONNNNVVVVEEEENNNNTTTTIIIIOOOONNNN
  180.        The matrices  are assumed  to be stored in a  oooonnnneeee ddddiiiimmmmeeeennnnssssiiiioooonnnnaaaallll CCCC aaaarrrrrrrraaaayyyy
  181.        in an analogous fashion as a Fortran array (column major). Therefore,
  182.        the element  A(i+1,j)  of matrix A  is stored  immediately  after the
  183.        element  A(i,j), while  A(i,j+1) is lda  elements apart from  A(i,j).
  184.        The element A(i,j) of the matrix can be accessed directly by reference
  185.        to  a[ (j-1)*lda + (i-1) ].
  186.  
  187.  
  188. AAAAUUUUTTTTHHHHOOOORRRRSSSS
  189.           Jack Dongarra, Argonne National Laboratory.
  190.           Iain Duff, AERE Harwell.
  191.           Jeremy Du Croz, Numerical Algorithms Group Ltd.
  192.  
  193.  
  194.  
  195.                                                                         PPPPaaaaggggeeee 3333
  196.  
  197.  
  198.  
  199.  
  200.  
  201.  
  202. ____SSSSYYYYRRRR,,,,____HHHHEEEERRRR((((3333FFFF))))                                                    ____SSSSYYYYRRRR,,,,____HHHHEEEERRRR((((3333FFFF))))
  203.  
  204.  
  205.  
  206.           Sven Hammarling, Numerical Algorithms Group Ltd.
  207.  
  208.  
  209.  
  210.  
  211.  
  212.  
  213.  
  214.  
  215.  
  216.  
  217.  
  218.  
  219.  
  220.  
  221.  
  222.  
  223.  
  224.  
  225.  
  226.  
  227.  
  228.  
  229.  
  230.  
  231.  
  232.  
  233.  
  234.  
  235.  
  236.  
  237.  
  238.  
  239.  
  240.  
  241.  
  242.  
  243.  
  244.  
  245.  
  246.  
  247.  
  248.  
  249.  
  250.  
  251.  
  252.  
  253.  
  254.  
  255.  
  256.  
  257.  
  258.  
  259.  
  260.  
  261.                                                                         PPPPaaaaggggeeee 4444
  262.  
  263.  
  264.  
  265.